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Abstract 

We calculate the diurnal and annual modulation of the signals in axion and 
WIMP dark matter detectors on Earth caused by a cold flow of dark matter 
in the Solar neighborhood. The effects of the Sun's and the Earth's gravity, 
and of the orbital and rotational motions of the Earth are included. A cold 
flow on Earth produces a peak in the spectrum of microwave photons in 
cavity detectors of dark matter axions, and a plateau in the nuclear recoil 
energy spectrum in WIMP detectors. Formulas are given for the positions 
and heights of these peaks and plateaux as a function of time in the course 
of day and year, including all corrections down to the 0.1% level of precision. 
The results can be applied to an arbitrary dark matter velocity distribution 
f(v) by integrating the one- flow results over velocities. We apply them to the 
set of flows predicted by the caustic ring model of the Galactic halo. The 
caustic ring model predicts the dark matter flux on Earth to be largest in 
December db one month. Nonetheless, because of the role of energy thresholds, 
the model is consistent with the annual modulation results published by the 
DAMA collaboration provided the WIMP mass is larger than approximately 
100 GeV. 

PACS number: 95.35.+d 
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I. INTRODUCTION 



The existence of cold dark matter (CDM), a species of as yet unidentified non-baryonic 
particle that makes up a large fraction of the energy density of the universe, was dramatically 
reaffirmed recently by the release [1] of the first-year data from the Wilkinson Microwave 
Anisotropy Probe (WMAP). WMAP confirmed the earlier discovery [2] of the accelerated 
expansion of the universe, and determined all cosmological parameters with unprecedented 
precision. The universe is found to be flat, with overall composition: 73% dark energy, 
23% dark matter, and 4.4% ordinary baryonic matter. [By definition, dark energy has large 
negative pressure (p/p ~ —1) whereas dark matter has negligible pressure (p << p).] The 
existence of dark matter had been inferred earlier from the dynamics of galaxies and galaxy 
clusters [3,4]. Dark matter is thought to be collisionless, and it must be cold [5]. The first 
attribute means that the dynamical evolution of the dark matter is the result of purely 
gravitational forces. We believe dark matter has this property because it reveals itself so 
far only through its gravitational effects and because it forms halos around galaxies. If dark 
matter were collisionful, it would dissipate its kinetic energy after falling onto galaxies and 
settle deep into the galactic gravitational wells, as ordinary matter does. The second at- 
tribute means that the dark matter must have sufficiently low primordial velocity dispersion 
to allow it to cluster easily on galactic scales. In contrast, neutrinos, although collisionless, 
have too large a velocity dispersion to be the dark matter in galactic halos [6] . 

Elementary particle physics has provided two well-motivated candidates for cold dark 
matter: the axion and the neutralino. Both were postulated for reasons wholly independent 
of the dark matter problem. The axion was postulated [7] to solve the strong CP problem of 
QCD, i.e. to explain why the strong interactions conserve P and CP in spite of the fact that 
the Standard Model as a whole violates those symmetries. It has the required cosmological 
energy density to be the dark matter if its mass m a is of order 10 peV [8]. In this mass 
range, it is extremely weakly interacting. The primordial velocity dispersion of dark matter 
axions is very small: 



where t is the age of the universe, and t the present age. The neutralino is a dark matter 
candidate [9] suggested by supersymmetric extensions of the Standard Model. Supersymme- 
try has been postulated to solve the hierarchy problem, i.e. to explain why the electroweak 
scale is small compared to other known scales such as the Planck mass, or the grand uni- 
fication scale. Neutralinos are examples of a broader class of cold dark matter candidates 
called WIMPs [10]. WIMP stands for 'weakly interacting massive particle'. The neutralino 
mass is of order 100 GeV. The primordial velocity dispersion of WIMPs is 



where mw is the WIMP mass. 

Fortunately, both dark matter candidates can be detected on Earth, at least in principle. 
As was mentioned already, CDM falls into the gravitational wells of galaxies and forms 
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invisible halos around their luminous parts. The evolution and structure of galactic halos 
can be modeled. Although there is at present no consensus on the validity of specific models, 
it is generally agreed that the density of CDM on Earth is of order 10 -24 gr/cm 3 . Dark matter 
axions on Earth can be detected [11,12] by stimulating their conversion to microwave photons 
in an electromagnetic cavity permeated by a large static magnetic field. The energy of the 
outgoing microwave photon equals that of the incoming axion. The spectrum of microwave 
photons in the cavity can be measured very accurately. So, if a signal is found, it becomes 
possible to measure the kinetic energy spectrum of CDM on Earth with great precision. 
Note that all CDM candidates have the same velocity distribution at any physical location 
because the velocities of cold collisionless particles are the outcome of purely gravitational 
interactions. WIMP particles can be detected on Earth by measuring the recoil energy of 
nuclei off of which WIMPs have scattered elastically in the laboratory [13]. Unfortunately, 
because the scattering angle is not measured, one cannot deduce the kinetic energy of the 
incoming WIMP on an event by event basis. 

Any signal of CDM in the laboratory will depend critically upon the details of the 
local halo velocity distribution [14]. The oldest and simplest model of galactic halos is 
the isothermal model [15]. It assumes that the dark matter particles form a self-gravitating 
sphere in thermal equilibrium. The model depends on only two parameters: the halo velocity 
dispersion and a core radius. It predicts that the local velocity distribution is Maxwellian 
in the non-rotating restframe of the Galaxy with velocity dispersion 5i>h a io — 270 km/s. In 
that case the flux of CDM particles on Earth is largest near June 2 [16], when the orbital 
motion of the Earth is most nearly in the same direction as the motion of the Sun around 
the Galaxy, and smallest six months later. 

However, the underlying assumption of the isothermal model, that present galactic halos 
are thermalized, is not justified. Even if a halo was thermalized during an early era of " violent 
relaxation" [17], the continuing infall of CDM particles into the galactic gravitational well 
soon adds a non-thermal component to the halo [18]. An isolated galaxy, such as our own, 
keeps growing by the infall of the surrounding dark matter. The dark matter flows which 
result from this continual infall do not thermalize over the age of the universe if gravity is 
the only force acting on the dark matter. It was shown in refs. [18,19] that the fraction of 
the local dark matter density in these discrete flows is large (of order one). The flows form 
peaks in velocity space [18] and caustics [21-23] in physical space. 

A^-body simulations [24] correctly describe the evolution and structure of galactic halos 
provided is large enough. How large must N be? The number of particles per galactic 
halo is of order 10 84 in the case of axions, and 10 68 in the case of WIMPs. Because they 
have negligible primordial velocity dispersion, CDM particles lie on a 3-dim. sheet in 6- 
dim. phase-space. The problem of understanding the evolution and structure of galactic 
halos is identical to the problem of understanding the folding of this phase-space sheet. 
The minimum number of folds of the phase-space sheet at the Sun's location in the Galaxy 
(equivalently, the minimum number of CDM flows on Earth) is of order 100 [18]. Because 
phase-space is 6-dim. it takes at least 100 6 = 10 12 particles to describe the phase-space 
structure of the Milky Way halo down to our location in it, 8.5 kpc from the Galactic 
center. This is a strict kinematical requirement. Present simulations have only of order 10 7 
particles and hence can at best give a crude description of galactic halos. Present simulations 
are also afflicted by two-body relaxation [25]: the simulated particles (each representing 
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approximately 10 7 solar masses) make hard two-body collisions with one another. (Two- 
body interactions between axions or WIMPs are of course entirely negligible.) Note that 
the formation of discrete CDM flows is evident in simulations [26] when appropriate care is 
taken. 

A model of the structure of isolated galactic halos [19,21,22,27] has been developed 
which explicitly realizes the fact that all CDM particles lie on a 3-dim. sheet in phase-space. 
Called the 'caustic ring halo model', it is described in some detail in Appendix B. It predicts 
a discrete set of flows at any physical location. The number of flows is location dependent. 
Caustic surfaces separate the regions with a differing number of flows. At the caustics, the 
dark matter density is very large. Evidence for the caustics predicted by the model has been 
found in observations of the Milky Way [28] and of other galaxies [29]. From fitting the 
model to the observations, a set of CDM flows at the Earth's location has been determined. 
The flows are listed in Table I. It is an update of a similar table published in ref. [27]. 
Because of our proximity to a caustic, a pair of flows dominates the local CDM distribution. 
It is given by the n = 5 entry in Table I. One member of the pair, dubbed the 'big flow', 
has of order three times as much density as all the other flows combined. 

It is not the purpose of our paper to review galactic halo models, or to try to settle the 
associated controversies. We described the main approaches to galactic halo modeling only 
to provide a context for our work, whose aim is much more modest. We wish to answer the 
question: given a single flow of dark matter of uniform density d and velocity v far from 
the Sun, what signals are produced in a dark matter detector on Earth as a function of time 
of day and of year? Our answer below takes account of the orbital and rotational motions of 
the Earth, and the gravitational fields of the Sun and the Earth. By including these effects, 
an accuracy of 0.1% in the answer is achieved. 

The effect of solar gravity [30] is such that a single uniform flow incident upon the Sun 
produces three flows downstream of the Sun, whereas there is only one flow upstream. See 
Fig. 1. The surface separating the one- flow region from the three-flow region is the location 
of a caustic, called the 'skirt'. It is a cone centered on the Sun, with opening angle equal to 
the maximum scattering angle of the particles in the flow. There is another caustic, called 
the 'spike', on a line downstream of the Sun, where the incoming flow gets focused by the 
gravity of the Sun. For the sake of clarity, we call the initial flow, with uniform density 
do and velocity Vq far from the Sun, the 'mother flow'. The mother flow produces near 
the Sun one or three position-dependent 'daughter flows' with densities di(r) and velocities 
Vi(r), where % — 1 or % — 1, 2, 3. The present axion and WIMP dark matter detectors are 
direction-independent. The signals produced in these detectors depend only on the densities 
di(t) and speeds v i(t) of the daughter flows in the laboratory rest frame as functions of time 
of day and of year. For given do and vq, we calculate di(t) and Vi(t) with the precision 
advertised above. 

We have been motivated by the expectation that the local CDM velocity distribution is 
dominated by a set of discrete flows, as in the caustic ring halo model. Our results can be 
applied directly to such flows. Our results can also be applied directly to local streams of 
dark matter caused by the tidal disruption of small satellite galaxies captured by the Milky 
Way [31]. Observations and modeling of the tidal disruption of the Sagittarius A dwarf 
galaxy by the gravitational field of the Milky Way indicates that the Sun is very close to 
that stream, and may be in it [32]. Finally, we stress that our results can be applied to any 
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local CDM velocity distribution f(vo) by simply integrating the single mother flow results 
over the initial velocity Vq. 

It was already mentioned that the isothermal model predicts the flux of dark matter on 
Earth to be largest near June 2, when the orbital motion of the Earth is most nearly in the 
same direction as the motion of the Sun around the Galactic center. The reason is that, in 
the isothermal model, the dark matter halo is not rotating. The dark matter particles in the 
isothermal model carry no net angular momentum. This is unrealistic because the angular 
momentum exhibited by the luminous matter is presumably caused by gravitational forces 
which must have also applied a net torque to the dark matter. The reasonable expectation 
is that the dark halo is rotating in the same direction as the visible matter. When halo 
rotation is included, it is not at all clear that the flux of CDM on Earth is largest on June 
2, and lowest six months later, as predicted by the isothermal model. In fact the caustic 
ring halo model, which includes net angular momentum of the dark matter particles in 
a systematic way, predicts that the CDM flux is largest in December ± one month. In 
the caustic ring halo model, the CDM flows of highest local density have velocities in the 
direction of galactic rotation and with magnitude (~ 480 km/s) larger than the speed (220 
km/s) of the visible matter [27]. In that case the annual modulation of the CDM flux on 
Earth has phase opposite to the prediction of the isothermal model. 

Motivated by this remark, J. Vergados [33], A. Green [34], and G. Gelmini and P. Gondolo 
[35] investigated the annual modulation in the caustic ring halo model of the signal searched 
for by WIMP dark matter experiments. In these pioneering papers, it was established 
that a flow of dark matter on Earth produces a plateau in the spectrum of recoil energies 
E r of nuclei off of which WIMPs have scattered elastically. The edge of the plateau is the 
maximum nuclear recoil energy E max for given WIMP mass m x and flow velocity Vp through 
the laboratory. E max is proportional to vp. The height of the plateau is proportional to Vp 1 . 
The total signal, obtained by integrating the recoil energy spectrum over recoil energies, 
is therefore proportional to vp ■ Vp 1 = vf, i-e. it is proportional to the flux of WIMPs 
through the laboratory as expected. So the total WIMP signal in the caustic ring halo model 
is largest in December, when the CDM flux is largest. However, as the above-mentioned 
authors [33-35] already emphasized, actual WIMP search experiments are generally only 
sensitive to a finite range of possible recoil energies, say from E\ to E 2 . Depending on the 
relative values of E\, E 2 and E max , the measured WIMP rate may be proportional to vf, or 
inversely proportional to vf, or something else altogether. 

The DAMA WIMP dark matter experiment in Gran Sasso Laboratory observes an annual 
modulation of its event rate, and interprets the effect as due to WIMPs [36]. The DAMA 
event rate is maximum on May 21 ± 22 days, in agreement with the modulation of the WIMP 
flux in the isothermal model (maximum on June 2). The region of parameter space claimed 
by DAMA appears inconsistent with the null observations of the CDMS [37], Edelweiss 
[38] and Zeplin [39] experiments. However, because the latter do not use the same target 
nuclei as DAMA, the comparison necessarily involves some assumptions about the WIMP 
couplings as well as assumptions about the structure of the Milky Way halo. We emphasized 
already the lack of consensus on the validity of any specific halo model. We also mentioned 
that the caustic ring model predicts that the CDM flux is largest in December ± one month. 
Such a prediction may at first sight appear inconsistent with the DAMA results but, because 
DAMA observes recoil energies over a finite range only, depending on the WIMP mass, and 
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for the reasons mentioned at the end of the previous paragraph, there is, in fact, no such 
inconsistency. A. Green concluded [34] that the caustic ring model, as described in [27], and 
the DAM A results are consistent if the WIMP mass is larger than 30 GeV. 

The present paper is guided by the following considerations. First, there is strong need 
for a detailed analysis of the annual and diurnal modulations of the signals in the axion dark 
matter experiment, with a precision commensurate with that of the experimental apparatus. 
If a signal is found in the ADMX cavity detector [12] presently taking data at Lawrence 
Livermore National Laboratory, the CDM kinetic energy spectrum will be measured readily 
with a precision of 10~ 5 , and this would soon be improved upon. It is desirable therefore to 
calculate the diurnal and annual modulation of the CDM density and velocity distribution 
on Earth with comparable precision. To this end we take account not only of the orbital 
motion of the Earth, but also the effect of solar gravity, the rotation of the Earth, and the 
effect of Earth's gravity. We use the previous results by two of us [30] on the effect of solar 
gravity. Second, since the axion and WIMP velocity distributions are the same, we present 
the results in a way which can be readily used by the WIMP dark matter experiments 
as well. The precision and sensitivity of the WIMP experiments is rapidly improving. If 
a candidate WIMP signal is found, the observation of the diurnal modulation and/or the 
effect of Solar gravity may provide an unambiguous confirmation. Third, we update the 
predictions of the caustic ring model in refs. [33-35] because the caustic ring model has 
evolved since those papers appeared. In ref. [28] the model was fitted to measurements of 
the Milky Way rotation curve and to a triangular feature in the IRAS map of the galactic 
plane. The IRAS feature is interpreted to be the imprint of the nearest caustic ring of dark 
matter on gas and dust in the galactic plane. If so, the Sun is very close to a cusp in the 
caustic ring nearest to us. This implies that the densities of the flows involved in that caustic 
(the n = 5 flows) are very large at our location. The fit requires an update of the list of 
local CDM flows predicted by the model. The new list of flows is given in Table I. The main 
difference with ref. [27] is that the fifth flows (n = 5) are now much more prominent. In fact, 
the revised caustic ring halo model predicts that the local dark matter velocity distribution 
is dominated by a single flow, the 'big flow', containing of order 75% or more of the local 
dark matter density. 

The remainder of this paper is organized as follows. In Section II we review the results 
of ref. [30] on the effect of the Sun's gravity on a cold collisionless dark matter flow. Sec. Ill 
describes the diurnal and annual modulation of the signals in the axion dark matter cavity 
experiment. Sec. IV describes the modulation of the WIMP signal. The predictions of the 
caustic ring halo model are compared with the DAMA results. We give our concluding 
remarks in Sec. V. Appendices describe our coordinate systems and the local CDM velocity 
distribution predicted by the caustic ring model. 

II. EFFECT OF THE SUN'S GRAVITY 

This paper analyzes the diurnal and annual modulation of the signals in dark matter 
detectors located on Earth due to discrete cold flows of dark matter in the Solar neighbor- 
hood. The annual modulation is caused by the Earth's orbital motion and by the spatially 
varying effect of the Sun's gravity on the flow. We summarize in this section relevant results 
from ref. [30] describing the effect of the Sun's gravity. 
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We assume the flow to be spatially and temporally uniform far upstream of the Sun 
with density d and velocity Vq. We assume also that the velocity dispersion of the flow is 
negligibly small. The effect of the Sun's gravity is such that there are three flows (i = 1, 2, 3) 
downstream, whereas there is only one flow (i = 1) upstream. See Fig. 1. To describe the 
flows, it is convenient to adopt cylindrical coordinates (z, p, 0) such that vq = VoZ, with 
< < 7r and — oo < p < +oo (instead of the usual < < 2n and < p < +oo). The 
range of z is — oo < z < +oo as usual. The flows at a given location (z, p, 0) have different 
impact parameters bi(z,p). We allow the b ri to have either sign as well. For flow 1, the 
impact parameter b\ has the same sign as p, whereas the impact parameters 6 2 and b 3 of the 
other two flows have opposite sign from p. See Fig. 1 . Flow 1 exists both upstream and 
downstream of the Sun, whereas flows 2 and 3 exist only downstream. Each flow is scattered 
through an angle where — O max < © < +9 max and max is the maximum scattering angle 
among all trajectories. Flows 2 and 3 are present downstream at angles — O max < 9 < +0 max 
and end on a conical caustic called the 'skirt', with opening angle 6 max . In the limit of a 
point mass Sun, @ max — > 7r, flows 1 and 2 exist everywhere, and flow 3 disappears because 
g?3 — > everywhere in that limit. Flow 3 always goes through the Sun. In addition to the 
skirt, there is a caustic located on a line downstream of the Sun, called the 'spike'. It occurs 
because the gravitational field of the Sun focuses the (collisionless) flow of dark matter. 

As a function of position (z, p, 0) the impact parameters of flows 1 and 2 are 

where a = and r = Vp 2 + z 2 . Eq. (2.1) assumes that the mass distribution of the Sun 

is spherically symmetric and that the particles comprising the flows in question did not go 
through the Sun at any point in the past. The latter condition is satisfied at most locations 
in the solar neighborhood, but not everywhere. Note that the quadratic equation: 

b 2 -bp- a(r + z) = (2.2) 

is solved by b\ and b 2 . 

The properties of the flows which have previously gone through the Sun can be expressed 
most readily in terms of the scattering angle 6(6 S ). b s is the impact parameter of the 
trajectory of scattering angle 6. We take b s to be always positive, the usual convention. 
The relation between and b s depends for small values of b s on the mass distribution inside 
the Sun. Q(b s ) was derived in ref. [30] for a particular solar model. The relation is two to 
one: for every value of the scattering angle there are two values of the impact parameter. 
Let us call them 6 s2 (6) and b s3 (Q) : with b s3 (Q) < b s2 (<d); see Fig. 2. Sufficiently far from 
the Sun, i.e. for r >> R Q where R Q = 6.961 • 10 10 cm is the solar radius, the scattering 
angle 6 of the particles in a flow that passed through the Sun is approximately the polar 
angle 9 = tan _1 (^) at that location. For such flows {i = 2, 3) 

\bi(6)\ = b st (Q = 9) + O(^) . (2.3) 

Far from the Sun, such as on Earth, the flows that passed through the Sun, and to which 
Eq. (2.3) therefore applies, are flow 3 everywhere and flow 2 for 9 near max . Fig. 3 shows 
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the impact parameters bi(0) for a particular solar model and a particular value of the initial 
speed Vq. Note that at 9 = flow 1 becomes flow 2 and vice versa. 

We now describe the flows by giving their density and velocity vector as functions of 
position. For all flows (i = 1, 2, 3), we may write the following expression for the velocity 



v Qi (r) = ±rv c 



'i + - - ( b -y - §v£ 



(2.4) 



Eq. (2.4) follows from energy and angular momentum conservation and is valid everywhere 
outside the Sun. It is also valid inside the Sun provided we replace ?f by ~ 2 ^ r ^ , where V(r) 

is the Sun's gravitational potential. The +(— ) sign pertains to where the flow is outgoing 
(incoming), i.e. where r > (r < 0). For the flows that did not go through the Sun, Eq. 
(2.4) may be rewritten in other useful ways (i = 1,2): 



v ei (r) = zv (l + — ) - p— (bi - p) 

= v + (l =F VTTy) (f - z) 

The densities of the flows that did not go through the Sun are 

1 



where 



di(r) 



Y 



do 
4 



VT+Y + 



VT+Y 



2a 



±2 



4a 

J [r + Z) ~ rsin 2 (f) 



(2.5) 



(2.6) 



(2.7) 



The spike caustic is evident because Y diverges as when p — > with z = +r, implying 
di,2 — \f^^f i n that limit. The densities of the flows that did go through the Sun are 

(i = 2,3) 



di 



b si (9) 



;i + o ( ^)) 



r 2 | sin ^ 

The skirt caustic is evident because ^ goes to zero at the maximum scattering angle 6 



(2.8) 



III. AXION SIGNAL MODULATION 

In this section, we discuss the diurnal and annual modulation of a signal in a cavity 
detector of dark matter axions due to a single cold flow of axions incident upon the solar 
system. We assume the detector is placed on Earth. 

The detector is an electromagnetic cavity permeated by a strong static magnetic field 
[11]. In this cavity galactic halo axions convert to microwave photons. The microwave 
power is amplified and spectrum analyzed. For every conversion, the photon energy equals 
the total (rest mass + kinetic) energy of the incoming axion: 
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= rn a c 2 = 2 / 1 u 2 3 u 4 \ 

Jl- (2)2 V 2V 8V / V ; 

where m a is the axion mass and v is the axion speed in the laboratory frame. The conversion 
is resonantly enhanced if the photon frequency f — j- falls within the cavity bandwidth. 
Dark matter axions are expected to have a mass in the 10 -5 eV/c 2 = 27r(2. 4179895 GHz)4r 
range with large uncertainties. Galactic halo axions have velocities of order 10~ 3 c. Hence 
the spectrum of microwave photons from a — > 7 conversion in the cavity detector has width 
of order 10~ 6 z/ a above v a = ^j^-. 

The microwave power from a — > 7 conversion is proportional to the quality factor Q of 
the cavity. Because it must be permeated by a strong magnetic field, the cavity is made of 
normal (i.e. non-superconducting) material. This allows Q ~ 10 5 . In signal search mode, 
the cavity is tuned over a wide frequency range, in the hope of finding u a . At a given 
frequency, the detector is sensitive to dark matter axions over a frequency range equal to 
the cavity bandwidth f /Q. Hence, if for example one octave in frequency (from / to 2f) is 
to be covered in one year, the amount t of time that can be spent at each cavity frequency is 
of order 100 sec. This in turn sets the maximum frequency resolution of the cavity detector 
in search mode to be of order 10 mHz. When this highest frequency resolution is achieved 
- as is presently the case with the ADMX detector [12] at Lawrence Livermore National 
Laboratory - the energy resolution 

AE vAv Af 



E c 2 / 
is of order 10 -11 and hence the velocity resolution 



(3.2) 



c 2 

Av ~ HT 11 10~ 8 c = 3 m/s . (3.3) 

v 

Note that if a signal is discovered the detector will remain at the axion mass frequency for 
much longer times, and the resolution will be correspondingly increased. 

Consider now the signal from a single cold flow of dark matter axions incident upon the 
solar system. It produces a peak in the frequency spectrum of microwave photons from 
axion to photon conversion in the cavity detector. The size of the peak is proportional to 
the density of the flow on Earth. The frequency of the peak is given by Eq. (3.1). The 
width of the peak is 

* / = / ^0.3.10-»J4) (3.4) 

where 5v is the velocity dispersion of the flow at the detector location. As the Earth orbits the 
Sun, both the density d and the flow velocity v change, producing an annual modulation of 
the signal. In addition, the rotation of the Earth produces a diurnal modulation of the signal 
frequency. We will discuss first the frequency modulation and then the size modulation. 

A. Frequency modulation 

We saw in Section II, that a single flow, of uniform density d Q and velocity v far from 
the Sun, produces one or three flows near the Sun depending on location, and hence one or 
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three peaks in the cavity detector of dark matter axions depending on time of year and on 
the direction of Vq relative to the Earth's orbit. Eq. (2.4) gives the flow velocities v Q i in the 
rest frame of the Sun. The velocities in the detector rest frame are 

Vi = v Qi - v orh (t) - v TOt (t) + Av i(B (3.5) 

where v orlJ (t) is the orbital velocity of the Earth, v rot (t) is the velocity component of the 
laboratory caused by Earth's rotation, and Av i(B is the effect of the Earth's gravity. We 
neglect here the effect of the gravitational fields of the other planets. The frequency shift 
due to Jupiter is of order AL ~ ~ 1CT 12 where Mj is Jupiter's mass and Rj its average 
distance from the Sun. Such effects are certainly measurable in the cavity detector but they 
are small and complicated to describe in detail. So we ignore them in this paper. The 
relativistic corrections in Eq. (3.1) are also of order 10~ 12 , and we neglect them as well. 

1. Effect of the Earth's gravity 

We have therefore 

Afi fi~V a -> -> A -> \2 

v a v a 2c 2 

= 7^( V ®i ~ ~ ^ot) + (3.6) 

where M e and i? e are the mass and radius of the Earth. The last equality in Eq. (3.6) follows 
from energy conservation. The effect of Earth's gravity is to increase all peak frequencies 
by the same amount 



Afi 



Va 



^ = 6.96 1(T 10 (3.7) 



independently of time. 



2. Effect of the Earth 's rotation 



To discuss the diurnal frequency modulation it is convenient to rewrite the first term on 
the RHS of Eq. (3.6) as 

— (v ei - Vrot) 2 = ^( V li + v rot ~ 2^ • £ot) ( 3 - 8 ) 

where 

%i = %i ~ ^orb (3.9) 

is the flow velocity in a reference frame which is co-moving but not co-rotating with the 
Earth, and 
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km 

v T ot = 0.465 — (— sin0 s x® + cos(p s y e ) cos/ 
s 



(3.10) 



is the velocity of the laboratory in that reference frame. We use Earth based coordinates 
(x®, y®, z®) such that i® is the Earth's rotation axis pointing North, x® points towards the 
Vernal Equinox (the direction of the Sun at the time of the vernal equinox) and y® = i® x x® . 
/ is the latitude of the laboratory 6 S = 2ttj- where t is sidereal time - i.e. the time elapsed 
since the Vernal Equinox last crossed the local meridian - and t s is the sidereal day. Eq. 
(3.10) neglects the asphericity of the Earth. The second term in Eq. (3.8) increases all peak 
frequencies by a tiny time-independent amount 



v. 



e rot 



rot 

2c 2 



^ = 1.2 10" 12 cos"/ 



(3.11) 



The last term in Eq. (3.8) can be spelled out: 

\ km 

— -v ei ■ v rot = 0.465 — cos/ (sin0 s v ei 



COS (J) s 



(3.12) 



Since t>®j is of order 300 km/s, the daily frequency modulation amplitude is of order 10~ 9 . 
Eq. (3.12) shows that the measurement of the amplitude and the phase of the daily frequency 
modulation determines the values of v®i X and v^ y on any day of the year. Moreover, the 
magnitude of v^ z is determined by measuring the day-averaged frequency shift away from 



day 
averaged 



~~ 2c 2 ^ V ® ix V ® iy v ® iz ^ ^ 



(3.13) 



The ambiguity in the sign of v®j Z is removed by observing the annual modulation; see below. 
Thus, on any given day, all three components of -u®j can be measured. 



3. Effect of the Sun's gravity and of the Earth's orbital motion 

Next we consider the annual frequency modulation, i.e. the time dependence of 
A/, 



v„ 



day 
averaged 



GM e 1 _ 2 

^--^fe-^orb) 



(3.14) 



Using energy conservation for the dark matter particles in the flow we have 



A/, 



GMr. 



1 _ 



c 2 R e 2c 2 c 2 r(t) 2c 2 

averaged 



orb 



(3.15) 



where t now stands for time of year. Furthermore, we may use energy conservation in the 
orbital motion of the Earth 



I 2 m _ ^0 _ GM e 
2 V ° a[ ) r(t) ~ 2a® 



(3.16) 



11 



to rewrite Eq. (3.15) as 



v„ 



day 
averaged 



2c< 



+ 



GMr,, 



r(t) 




1 - 



(3.17) 



Here, a e = 1A.U. = 1.4960 10 13 cm is the semi-major axis of the Earth's orbit. 

Let (r, 0, $) be spherical coordinates centered on the Sun. The polar angle is measured 
relative to Z, the direction perpendicular to the plane of the Earth's orbit, and $ = in the 
direction of the Vernal Equinox, i.e. the direction of the Sun as viewed from Earth at the 
time of the vernal equinox. Thus the Earth is at $ = at the time of the autumnal equinox 
(near Sept. 23). The angular position $(t) of the Earth and its distance r(t) to the Sun are 
given by the parametric equations: 



tan( ^ ) = Vi^; tan 2 

r = a e (l - e© cos^) 



t t Ti 



P \GM & 



(ip - e© simp) 



(3.18) 



where e© = 0.0167 is the eccentricity of the Earth's orbit and ip is the 'eccentric anomaly' 
parameter. In Eqs.(3.18) the subscript p indicates that the quantity is evaluated at peri- 
helion, which occurs near Jan. 4. Using Eqs. (3.18), the second term on the RHS of Eq. 
(3.17) may be rewritten as 



GMq 



+ 2e ffi cos 



with 



GM G 



© 



c 2 a a 



= 0.987 10" 



(3.19) 



(3.20) 



Eq. (3.19) describes an increase of all frequencies by the same amount, the time- dependent 
part of which is of order 10~ 10 . 

Finally we analyze the last term in Eq. (3.17). Although we kept it for last, it is by far 
the largest contribution to the frequency modulation. For the sake of clarity, we first neglect 
the eccentricity e© of the Earth's orbit. 



4- Neglecting the eccentricity of the Earth's orbit 

= limit 

(2 = 1,2,3) 



In the e ffi = limit, v orh = Combining this with Eq. (2.4) we have for all flows 



-^v Qi -v orh = +6-$-J—2-b i , (3.21) 



<r <r v a'a 
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since $ • f — 0. Let (9 , $o) be the polar coordinates of the flow velocity v = v z in the 
Sun centered coordinate frame defined above. We have then (for p > 0) 



- - f cos 6* — z - £ • $ 

e-$ = — — — $ = — — 

sin sin 

sin6 sin($ - $ ) 



\/l - sin 2 6 cos 2 ($ - $ ) 
In terms of a = GA i & , we have therefore 



(3.22) 



1^ vlfabi sme sin $- $ ) ,„ OQ x 

- — v Qi ■ v orh = + — J , (3.23) 

c c V a © a © - sin 2 6 cos 2 ($ - $ ) 



where $ = 2-7T— — , and £ is time since the last autumnal equinox. 

year ' ^ 

Let us estimate the frequency modulation amplitudes predicted by Eq. (3.23). For flow 1, 
|&i| ~ a ffi generically, i.e. away from the spike caustic. Since 



a ^ = 1 .4746 10»cm , (3.24) 



the frequency modulation amplitude is typically of order 1CT 7 for flow 1. For flow 2 [see 
Eq. (2.1)], \b 2 \ ~ a ~ 10~ 2 a e generically, i.e. away from the spike. Hence the frequency 
modulation amplitude is typically of order 1CT 9 for flow 2. Close to the spike, flows 1 and 
2 are alike: bi ~ ±& sp ikc with 

/300km/s s 



&spikc = ^aa® = 2.1005 10 12 cm ^ ^T/ ' ( 3,2S ) 

Therefore, near the spike, the frequency shifts for flows 1 and 2 are both of order 10~ 8 . For 
flow 3, 63 ~ Rq everywhere, and its frequency modulation amplitude is therefore of order 
10- 9 . 

Let us consider more closely the frequency shifts of flows 1 and 2 when the Earth comes 
near the spike caustic, i.e. when Y = . 2 " g . >> 1. For the Earth to come near the spike, 

1 rsm^(|) 1 ' 

the unperturbed flow direction z must lie close to the ecliptic plane. So let 60 = \ — 5, and 
$ - $0 = e, with 5, e « 1. Then 

■ v OTh = ±V2 °^ & 6 , (3.26) 
c 2 U2 c 2 a e y/5 2 + e 2 

where we used Eqs. (3.23) and (3.25). 

5. Including the eccentricity of the Earth 's orbit 

Finally, let us include the corrections due to the finite eccentricity of the Earth's orbit. 
Using Eqs. (3.18), we have: 
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\GMr. 







a® 1 — e e cos^ 
Combining this with Eq. (2.4), we find (i = 1, 2, 3) 



e© sin ip f + Jl - e% $ 



1 „ „ 

~~^ v Qi ' "^orb 



c 2 r 



± e e sin^M/ 1 + 



2a 



-A 



.2 u * 



sin6 sin($ - $ c 



^/l - sin 2 O cos 2 ($ - $ ) 



(3.27) 



(3.28) 



In Eq. (3.28) the ± sign corresponds to flows with r < 0. Flows 2 and 3 always have f > 0. 
$, ^ and r are given as functions of time t by Eqs. (3.18). The impact parameters hi are 
given in Eqs. (2.1) and (2.3), with z = rcos9, p = rsin#, and 



cos 9 = sin ©o cos($ — $ ) 



(3.29) 



6. Frequency shifts on the measurement time scale 

Finally, we estimate the frequency shift of a peak in the axion signal in a 100 sec time 
interval. This is of order the measurement integration time of the cavity detector in search 
mode. The shift due to the Earth's orbital motion is at most of order 

2ix 

A/ 1 orb ~10" 7 / (100 sec) ~ 2 • 10" 12 / (3.30) 

ioo B ec year 

for flow 1, less for the others. The shift due to the Earth's rotation 

A/1 rot £ 1.5 • 10~ 9 — (100 sec) ~ 10" 11 / . (3.31) 
J 1 iooscc day v /j v i 

We see that the effect of Earth's rotation is larger than the effect of its orbital motion on 
this short time scale. In the ADMX experiment, / ~ 1 GHz and t ~ 100 sec [12], so 
that A/ ~ 10 mHz, and hence the frequency resolution Sf = j ~ 10 mHz is as high as is 



profitable in search mode. 



B. Density modulation 

The size of a peak in the axion signal is proportional to the density of the corresponding 
flow at the detector location. A peak's size is much less accurately measured than its 
frequency. So we neglect in this subsection the subleading effects associated with the Earth's 
gravity and rotation, and the eccentricity of the Earth's orbit. 

The density of the flows that did not go through the Sun is given by Eqs. (2.6) and (2.7), 
with r = a®, $ = 2n-^, and 6 given by Eq. (3.29). Away from the spike, i.e. for Y « 1, 
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d, • d d 2 d (^-Y 2 + 0(F 3 )) ~ ^ 1 " V 



16 4 \ a(B J sin 4 (|) 

.do 2.43 10- (^^V * • (3.32) 
\ v J sm 4 (§) 



The above approximation breaks down when Y ~ 1 or 



< _ o /300 km/s N 



e ~ i6° y — v^) ' ^ 3 ' 33 ^ 

Close to the spike (Y >> 1), the densities of flows 1 and 2 are large and nearly equal: 

do PT7 , /OOffi 1 , / a 1 ,„„,, 

4 V 2 p V a © ^2(52 + e 2 ) 

with <5 = | — 6 and e = $ — <3> . Note that we are assuming here that the spike is truly a line 
caustic. The spike is a line caustic if the incoming flow is uniform and the Sun is spherically 
symmetric. However, a line caustic is a degenerate case. Generically, dark matter caustics 
are surfaces. When observed with sufficient resolution, the central line of the spike will be 
seen to be a caustic surface whose properties depend on the asphericity of the Sun and the 
lack of uniformity of the incoming flow. 

The density of the flows that did go through the Sun are given by Eq. (2.8). Let us first 
consider the density of flow 3 downstream of the Sun. Near 6 = 0, b^{9) = —B9 where B 
depends on the mass distribution inside the Sun; see Fig. 2. Hence 

d 3 (6) ~ d (—) (3.35) 

for 6 « 1. For the model considered in ref. [30], B ~ 0.7R Q for t> = 300 km/s, and 
hence d 3 ~ 10~ 5 do. Generically, i.e. away from the spike, d 2 and d 3 have the same order of 
magnitude. 

Finally we consider the density of flows 2 and 3 near the skirt caustic. Let b c be the 
value of the impact parameter b s for which the scattering angle G is maximum. For G near 
G ma x 

0(6 B ) = G max - l -H{b s - b c f (3.36) 



with 

d 2 G 



H 



dbl 

Inserting this into Eq. (2.8) we have 



> . (3.37) 

b s =b z 



d 2 ~ d 3 ~ ^ (3.38) 

r 2 sinG max y2/d(G max -#) 
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for 9 < ma x- For 9 > 6 max , d 2 and d 3 vanish. The densities of flows 2 and 3 diverge for 
9 — > @ max ,- a s the inverse square root of distance to the caustic surface. This behavior 
is characteristic of a simple fold catastrophe. However, the coefficient of the singularity is 
inversely proportional to r 2 . At the Earth's distance to the Sun, the skirt caustic is quite 
feeble, as we will see. 

The values of b c , O m ax and H depend on the solar mass distribution and on vq. For the 
solar mass distribution considered in ref. [30] and t> = 300 km/s, one has b c ~ 1.35i? , 
©max — 107° and H ~ in which case 

1 53 10~ 6 

d 2 ~d 3 ~ d ' = (3.39) 

V "max — 9 

for r = a e . In view of Eq. (3.29), the Earth crosses a skirt caustic twice per year if 

| cos m ax | < sin ©o . (3.40) 

The skirt crossing times are t sk = "^sk^ff 1 where $ s k are the solutions of 

cos ©max = sin @ cos($ sk - $ ) ■ (3-41) 

If sin@o < | cos ©max |, the Earth never crosses the skirt caustic. In the latter case, flows 2 
and 3 are always absent (present) on Earth if cos© m ax > (<) 0. As a function of time 



d b c 1.58 10 3 / sec 

<M*) - 777 ■> / = = = = \ u 7 i ' ( 3 - 42 ) 



v^al ^sinQo sin© m ax | sin($ sk - $ )| V I* ~ 1 



near the time of caustic crossing, on the side with two extra flows. For some flows the 
orientation of the Earth's orbit may be such that |@ ± @ m ax| — tt/2 in which case the 
Earth's trajectory grazes the skirt caustic, and the density enhancement lasts a relatively 
long time. 

IV. WIMP SIGNAL MODULATION 

In this section, we review briefly the equations describing WIMP detection in the labora- 
tory, derive the signal modulation for a single flow with zero velocity dispersion, and discuss 
the implications for the caustic ring halo model. 

A. General formalism 

The experimental approach to WIMP detection in the laboratory is to measure the 
recoil energy E r of nuclei off of which WIMPs have scattered elastically [13]. See ref. [14] 
for reviews. Although it is highly desirable to measure the recoil direction as well, this 
capability has not yet been achieved in present detectors. 

The differential event rate per unit detector mass is 

dR=—^vf(v)-^ I d\q\ 2 dv (4.1) 
m x m N M ; d\qf m 1 ' 
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where m x is the WIMP mass, d the local WIMP mass density, mjy the target nucleus mass, 
f(v) the normalized CDM velocity distribution in the laboratory reference frame, q the 
momentum transfer, and da/d\q\ 2 the differential cross-section for elastic WIMP-nucleus 
scattering. The recoil energy deposited in the detector is 

■ < 4 - 2 > 

For given initial WIMP velocity v, the maximum momentum transfer (and hence maximum 
recoil energy) occurs when the scattering is back-to-back in the center of mass frame, in 
which case \q\ = |flL = = 2m r v, where m r = mjv , m ' x is the reduced mass. The differential 
cross-section is usually expressed in terms of a "standard" cross-section at zero-momentum 
transfer, called <7 , and a form factor F(\q\), normalized so that F(0) = 1: 

The model-dependent factors describing the interactions of a specific WIMP enter oo, while 
the nuclear effects are encoded in F(|(fj). 

Integrating over the WIMP velocity distribution, we obtain the event rate per unit recoil 
energy: 

dR cr d ^,,^2 f°° f( v 



dE r 2m x m 2 r 



(Tod ^n-A\2 T(E r 



y/7rm x m% 220 km/s 



where \q\ = \/2m^E r , and v m - m = is the minimum velocity needed to transfer a given 
energy E r . T(E r ) is a dimensionless function containing all relevant information about the 
velocity distribution. It is conventionally defined as 

T(E r ) = ^(220 km/s) f°° ^ dv . (4.5) 
2 Jv min v 

Finally, the total event rate within an observed range of recoil energies is 

r E i dR 

R= ^dE r (4.6) 
J El dE r K J 

where E\ and E% are respectively the minimum and maximum values of that range. 



B. The single flow plateau 

Consider a single cold (i.e. zero velocity dispersion) flow, of density dp and speed vp 
relative to the laboratory. In that case 

f(v) = 6(v - v F ) (4.7) 
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and 



§r = 9 F\\V\) H(v F - v mm ) (4.8) 

dE r 2m x mjv F 

where H(x) is the Heaviside step function. In the limit where the form factor is slowly 
varying over the relevant momentum transfer range, the contribution of a single cold flow 
to the recoil energy spectrum is a flat plateau [33-35,31]: 



dR oo^i 



for E r < E max 



dE r 2m x mj vf 

= for E r >E max (4.9) 



where 



£ max = ^ (4.10) 
m N 

is the maximum recoil energy for the flow. The total rate within an observed range of recoil 
energies, E x < E r < E 2 , is 

R= n [^(mm{E 2 ,E max })-^(mm{E 1 ,E max })\ (4.11) 

2m x mf,VF 

where 

HE) = T F\\q\) dE r . (4.12) 
Jo 

For Ei = and E 2 = oo, we recover the usual expression for the total rate (per unit detector 
mass) 

R = (TtotVF-^— (4.13) 

where 

1 f E 



^tot = v 01 ^- I F (\V\) 2 dE r ( 4 - 14 ) 

-'-'max •'O 

is the total cross- sect ion. 

The differential rate ^=r- as a function of E r is a plateau whose height is inversely pro- 
portional to the flow velocity vp, and whose edge is proportional to v F . So the total rate 
increases with increasing flow velocity proportionately to vp, as one would expect, but the 
rate within an observed recoil energy range {E\ < E r < E 2 ) decreases with increasing flow 
velocity if E 2 < E mgx . 
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C. Signal modulation for a single flow 



As discussed in Section II, a single cold homogeneous flow of density d and velocity v 
far upstream of the Sun, produces in the vicinity of the Sun one or three flows, depending 
on location. Formulas for the densities di(r) and velocities v &i (f) of the flows (i — 1 or 
% = 1,2,3) were given in Section II for all f. The flow densities di(t) and speeds Vi(t) on 
Earth as functions of time of day and of year were obtained in Section III. (The latter section 
is motivated by the axion dark matter experiment. However the formulas can be applied 
to WIMPs as well.) To obtain the time- dependent signal -jf§-(t) caused by a single mother 
flow, one must therefore replace dp with di(t) and vp with Vi(t) in Eq (4.8), and sum over i. 

The expressions for di(t) and Vi(t) in Section III, however, include many corrections to 
which present WIMP detectors are not sensitive. In particular, flow 3 is unlikely to be 
observed in present WIMP detectors because its density d^ is of order 10 _5 <io- Likewise, 
flow 2 is unlikely to be observed except when the Earth approaches the spike caustic. The 
most important effect for WIMP detectors is the annual modulation of the velocity vi of 
flow 1 due to the Earth's orbital motion. It is of order 10%. A short list of the lesser effects 
is as follows. The Earth's gravity increases all the Vi by a constant amount of order 0.1%. 
The Earth's rotation produces a diurnal modulation of the Vi of order 0.3%. The effect of 
the Sun's gravity on the Vi is of order 1%. The effect of the Sun's gravity on the di is of 
order 10 -4 except near the spike caustic where di and d 2 diverge. The corrections due to 
the eccentricity of the Earth's orbit are of order 0.1%. 

Here, we recapitulate our results for di(t) and Vi(t) including all effects down to the 0.3% 
level of the diurnal modulation. So we neglect the eccentricity e e of the Earth's orbit. We 
also neglect the existence of the third flow and all other effects, such as the skirt caustic, 
associated with the finite size of the Sun. With these approximations, the densities are given 
by Eq. (2.6) with 

Y = a(—) 2 . (4.15) 

V v J 1 - sin6 cos($ - $ ) 



We used - = — = (i^-) 2 since the eccentricity of the Earth's orbit is neglected. As before 



(0o, $o) are the solar polar coordinates of the initial flow direction -Oq, and $ = 27r^: where 



t is time since the autumnal equinox (near Sept. 21). The flow speeds on Earth are 



Vorh-V rot ) 2 (4.16) 



where v & i = v + Av &i is the flow velocity corrected for the effect of solar gravity, and v or h 
and v rot are the velocity components the laboratory has in the rest frame of the Sun due to 
the Earth's orbital and rotational motions. We have in solar coordinates 

v rb = v OT h (— sin $ X + cos $ Y) 

v TO t = frot cos / [— sin <fi s X + cos s (cos 7 Y — sin 7 Z)] 

v Q i = \v Q {\ ± VTTY) [sin O (cos $o X + sin $ Y) + cos O Z] 

+ _ Uo (i T Vl + Y) (cos $ X + sin$f) , (4.17) 
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where t> or b = 29.8 km/s and i> ro t = 0.46 km/s. I and <p s were defined in Section III, after Eq. 
(3.10). 7 = 23.44° is the angle between the Earth's equatorial plane and its orbital plane. 
The latter is usually referred to as 'the ecliptic'. 

The most important effect is caused by the Earth's orbital motion. If we neglect the 
lesser effects of solar gravity and terrestrial rotation, we have d\ = do, d 2 = 0, and 



vi = y [y - v orh ) 2 ~ v - z- v orh 

^v + v orh sin6 sin($ - $ ) . (4.18) 

At this order, the speed modulation is sinusoidal with amplitude sin O = 

10% ( 3 °°^ m ^ s ) sin ©o- The modulation amplitude is largest when the initial flow direction 
is near the ecliptic plane. In that case, however, the spike caustic is near the ecliptic plane 
as well and the effect of solar gravity is important when the Earth approaches the spike 
caustic. 

If we include the effect of solar gravity, but still neglect that of Earth's rotation, we find 
(vi) 2 = v 2 + 3 Vo 2 rb - v Q -v orh {l ± VT+Y) , (4.19) 

where we used vl rh = Solar gravity increases all flow speeds by the relative amount 

( 3 ^ k ) 2 = 1% ( 3 °°^ m/s ) 2 an d modifies the z ■ v orh term in a position-dependent way. When 
approaching the spike, Y diverges. However, Vi remains finite because vo'V rb ~~ > there. On 
the other hand, as Eq. (2.6) shows, both d\ and d 2 diverge at the spike. Setting O = | — 5 
and $ = $ + e, we have 



dx ~ d 2 ~ do (4.20) 
for e, 5 « 0. The maximum value of the density enhancement F = j^{d\ + d 2 ) is 



~ V2 ^ i . (4.21) 
v 6 



The duration of a peak with enhancement factor larger than F is 



At= ^S J(%i)2-1 (4.22) 



7T 



for 1 << F < F max . 

Finally, including the effect of Earth's rotation, we have 



= (%i - #orb) 2 + 2GM& - 2 ({J oi - V orh ) ■ V Iot 

dm 



V + 3 V orb - ^0 • #orb(l ± V 1 + Y) - 2 v Qi ■ v rot + 2 v orh ■ v mt . (4.23) 
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The last term 

2 v orb • v TOt = 2 v OT bV ro t cos I [cos 7 cos $ cos (p s + sin $ sin0j (4.24) 

is a correction of order 1CT 3 . The penultimate term in Eq. (4.23) is 

-2v &i ■ v TOt = -v Q v rot cosl [(1 ± VT+Y)(- sin0 s sin© cos$ 
+ cos 7 cos S sin6o sin$ — sin 7 cos0 s cos O ) 
+ (1 =F Vl + Y)(— sin0 s cos$ + cos 7 cos0 s sin$)] (4.25) 

with Y given by Eq. (4.15). 

D. WIMP signal modulation in the Caustic Ring Halo Model 

The caustic ring halo model predicts a set of discrete flows in the solar neighborhood. 
The model is briefly described in Appendix B. The list of flows is given in Table I in galactic 
coordinates, and in Table II in solar coordinates. The local dark matter distribution is 
dominated by one pair of flows. They are associated with the caustic ring of dark matter 
nearest to us, and are given by the n = 5 entries of the Tables. One of the n = 5 flows, dubbed 
the 'big flow', has density of order three times the sum of all the other flows combined. This 
feature simplifies the analysis of the model's predictions for the WIMP signal modulation. 
One complicating factor is our present ignorance whether the velocity of the big flow is 
in which case the velocity of the other n = 5 flow is v£, or vice versa. In any case, neither 
v§ nor v£ lie close to the ecliptic plane, and hence the Earth never approaches the spike 
caustic of the big flow. 

Fig. 4 shows the speeds on Earth of the daughters of the big flow as a function of time 
of year if the big flow initial velocity is v§. The Earth enters the skirt near March 30 and 
exits it near Dec. 15. So flows 2 and 3 only exist between those two dates. They are at 
any rate very hard to observe in a WIMP detector since their densities are suppressed by 
five orders of magnitude relative to flow 1. Flow 1 has speed modulation of order 7%. The 
speed, and hence the flux on Earth, is largest near Nov. 10 and lowest near May 10. 

Fig. 5 shows the annual speed modulation of the big flow if its velocity is v$. In this 
case also, flow 1 has a speed modulation of order 7% but the maximum occurs near Jan. 20 
and the minimum near July 20. 

Table II shows two flows whose initial velocity vectors vq are close to the ecliptic plane, 
namely 4+ and 10+. For reasons mentioned in Appendix B, the uncertainties on the flow 
velocities v n± are larger for n < 4 than for n > 5. So, flow 4+ is less certain to be close to 
the ecliptic plane than flow 10+ and therefore less certain to show spiky behavior. Figs. 6 
and 7 show the velocity and density modulation of flow 10+. The behaviour near the spike 
is, of course, sensitive to the precise value of v^ + . As an example, we chose v l % + = 2.6 km/s 
in Figs. 6 and 7. The Earth approaches the spike near Dec. 20. The densities of flows 1 and 
2 have a large temporary enhancement (-F max ~ 16), and their velocities change abruptly 
then. 

The WIMP signal modulation for the caustic ring halo model is dominated by the speed 
modulation of the big flow. This is evident in Figs. 8 and 9, which show T(E r ) at those 
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times when the big flow speed is largest and smallest. In Fig. 8 the densities with 
n = 5,6,7,8,9 are all assigned the velocities v^, whereas in Fig. 9, they are all assigned 
the velocities v^. Each flow modulates as described in the previous subsection, with its 
individual phase and amplitude determined by its initial flow direction. 

To give a qualitative discussion let us neglect all but the big flow. To include the other 
flows is not difficult, but it complicates matters considerably without changing the results 
much. The annual modulation depends sharply on the range, E\ < E r < E 2 , of measured 
recoil energies. The edge of the big flow plateau as a function of time of year is 



£max *) = — (vf) 

m N 



„ 60 km/s . ^ , . / 2?r , 

1 + + ^ sin Of sin (t-t±) 

v 5 ' Wear 



(4.26) 



with either v$ = 255 km/s, 65 = 41° and t_ = t on Aug. 10, or vt = 265 km/s, &t = 36° 
and t + = t on Oct. 20, depending on whether the initial velocity of the big flow is ih~ or 
v£. The height of the big flow plateau modulates with ~ 7% amplitude whereas the edge 
modulates with ~ 15% amplitude. The plateau height and edge modulations have opposite 
phase. See Eqs. (4.8-4.10). 

If E 1 and E 2 are both less than the minimum value -E max ,min of E max (t) in the course 
of the year, the observed event rate modulates with amplitude of order 7%, and maximum 
near either May 10 or July 20. This prediction of the caustic ring halo model is independent 
of the nuclear form factors involved. Such a modulation is very similar to that which occurs 
in the isothermal halo model, which predicts 7% modulation with the maximum occurring 
near June 2. This is so in spite of the fact that the isothermal and caustic ring models could 
hardly be more dissimilar. 

For other values of E\ and E 2 , however, the rate modulation in the caustic ring halo 
model is very different. If E\ and E 2 are both larger than the maximum £ max ,max of E max (t) 
in the course of year, the WIMP rate vanishes. If E\ is less than -E max ,mm and E 2 is larger 
than £max,max, the annual modulation amplitude is 7% (if E\ = 0) or larger, with the 
maximum occurring near November 10 or January 20. If E 1 and/or E 2 fall between E maxmin 
and -E'max.max, the time derivative of the WIMP event rate changes abruptly when E max (t) 
equals E\ or E 2 . In particular, if -E max , m in < E\ < E max>max , the signal is on part of the year 
and off for the remainder. The different kinds of behavior are easy to derive by combining 
Eqs. (4.9), (4.10) and (4.26). 

The most dramatic prediction of the caustic ring model is that, for certain values of 
the recoil energy E r , those between -E max , m in and -E maXimax , 4§- from the big flow (of order 
75% of the total local density) is zero for a definite period of the year, and non-zero for the 
remainder. 

A question of particular interest is whether the caustic ring halo model is consistent with 
the annual modulation observed by the DAMA experiment [36] . DAMA searches for WIMP 
dark matter by measuring the scintillation produced in Nal crystals when WIMPs scatter 
off nuclei of the crystal. Of order 30% ("quenching factor" q^ a = 0.3) of the nuclear recoil 
energy is measured as scintillation if the WIMP scatters off a Na nucleus, whereas of order 
9% is measured if the WIMP scatters off a I nucleus (qi = 0.09). Although Na has the larger 
quenching factor, in most particle physics models the WIMP elastic scattering rate in Nal 
is dominated by scattering off I nuclei because the latter are heavier (Aj = 127 vs. A Na = 
23). 
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DAMA observes an annual modulation of the scintillation rate with amplitude of order 
0.02 counts/day kg keV ee in the energy range 2-6 keV ee . The ee subscript, which means 
'electron equivalent', indicates that the scintillation energy is being given, not the nuclear 
recoil energy. The modulation is consistent with being sinusoidal, with phase such that the 
maximum rate occurs on May 21 ± 22 days. No significant annual modulation is observed 
in the 6-14 keV ee energy region. 

The annual modulation observed by DAMA in the 2-6 keV ee region is consistent with 
the caustic ring halo model provided E max is large enough. If the WIMP-nucleus interaction 
is spin-independent, the signal is dominated by scattering off I nuclei. The requirement that 
Qi E max > 6 keV implies that the WIMP mass is larger than about 190 GeV. WIMP masses 
as small as 100 GeV may be accommodated depending on the details of cross-sections, 
quenching factors, and nuclear form factors. Fig. 10 shows the event rate modulations in 
the range 2-6 keV ee for the caustic ring halo model with the big flow velocity equal to v§ 
for m x = 120 GeV and m x = 50 GeV, assuming the WIMP-nucleus interaction is spin- 
independent and the nuclear form factors are equal to one. 

The lack of annual modulation observed by DAMA in the 6-14 keV ee range, given that 
DAMA does observe a modulation in the 2-6 keV ee range, may be hard to accommodate 
in the caustic ring model. It would require some degree of cancellation among competing 
effects, e.g. a modulation with one phase during part of the year and with the opposite phase 
during the remainder. This can occur if qj E mSLX is close to 14 keV. However to compare 
the model with the DAMA data in such scenarios one would have to know the details of 
the nuclear form factors and make some assumptions about the particle physics. Such a 
comparison is beyond the scope of this paper. 

Instead let us emphasize that DAMA can test the caustic ring halo model by measuring 
-j^-(E r ,t) and asking whether for some values of E r , of order 75% of the event rate is off 
part of the year and on for the remainder, as predicted by Eq. (4.26). 

V. CONCLUSIONS 

For an arbitrary single cold flow, dubbed a 'mother flow', of initial velocity vq and density 
d far from the Sun, we derived the densities di(t) and speeds Vi(t) on Earth of the daughter 
flows as a function of time of day and time of year. We included all corrections down to the 
10~ 3 level of accuracy. In order of decreasing importance, the following effects are relevant: 
the orbital motion of the Earth (10%), the effect of solar gravity (typically 1% but of order 
100% or even larger if the Earth approaches a caustic), the Earth's rotation (0.3%), and the 
Earth's gravity (0.1%). Although we did not systematically study the next set of corrections, 
it seems that they are all much smaller. We expect the asphericity of the Earth and that 
of the Sun to produce corrections of order 10" 5 , and we already mentioned that Jupiter 
introduces corrections of order 10~ 6 . So, it may be that our calculations are in fact accurate 
at the 10 -5 level of precision. 

Our results are directly applicable to discrete dark matter flows and streams. However, 
they can also be applied to an arbitrary, continuous and/or discrete, dark matter velocity 
distribution f(v ), simply by integrating over v . 

The flows predicted by the caustic ring halo model were discussed as an example. They 
are listed in Tables I and II. The model predicts that the local velocity distribution is 
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dominated by a single flow, dubbed the 'big flow'. It is one of the flows associated with the 
caustic ring of dark matter nearest to us. There is a two-fold ambiguity in the velocity of 
the big flow. In one case (v 5 '), the big flow initial velocity has magnitude 255 km/s, and 
is pointing in the direction of declination 5$ = 31° and right ascension = 20. 2 h . In the 
second case (■w 5+ ), its initial velocity has magnitude 265 km/s, and direction 5$ = 60°, a£ = 
23.8 h . 

In the cavity detector of dark matter axions a cold flow produces a narrow peak in the 
spectrum of microwave photons from axion conversion. A peak's frequency can readily be 
measured with an accuracy of 1CT 11 . This implies an accuracy of order 10~ 5 in the flow 
speeds. [There is no difficulty in measuring a peak's frequency with 10~ 14 accuracy using 
atomic clocks. However, in the 100 sec signal integration time typical of the experiment in 
search mode, a peak shifts by an amount of order 10~ n ; see Eq. (3.31)]. The existence 
of the big flow allows the axion experiment to be made more sensitive by searching for the 
corresponding narrow peak. This capability has already been implemented in the ADMX 
experiment [12]. If a signal is found, it will be possible to measure for each mother flow the 
densities and speeds of the three daughter flows year round. Since the properties of flows 2 
and 3 near the skirt caustic depend on the mass distribution inside the Sun, the latter may 
be investigated. 

In a WIMP detector, a cold flow produces a plateau in the recoil energy spectrum 
of nuclei off of which WIMPs have scattered elastically. The edge of this plateau is the 
maximum recoil energy for given flow speed [see Eq. (4.10)]; it is proportional to the flow 
speed squared. On the other hand, the height of the plateau is proportional to the (flow 
speed) -1 . The annual modulation due to a single cold WIMP flow depends on the range 
of observed recoil energies. If the whole observed range is below the plateau edge, the 
event rate modulation has the opposite phase and the same relative amplitude as the speed 
modulation. If the whole plateau is in the observed range, the event rate modulation has 
the same phase and the same relative amplitude as the speed modulation. If only part of the 
plateau is in the observed range, the behaviour is more complicated. The speed modulation 
of the big flow has amplitude of order 7%. The maximum speed occurs near May 20 if the 
big flow velocity is -u 5- , July 20 if the big flow velocity is v 5+ . 

The most dramatic prediction of the caustic ring model for WIMP searches is that, for 
a certain range of recoil energies, the event rate due to the big flow (of order 75% of the 
signal) is on during a predicted part of the year and off for the remainder. 

The DAMA results published so far are not inconsistent with the caustic ring halo model. 
The modulation observed in the 2-6 keV ee range is consistent with the model provided the 
WIMP mass is large enough (~ 100 GeV). The lack of observed modulation in the 6-14 
keV ee range may be inconsistent with the model, but to actually show this one needs to know 
the nuclear form factors and to make some assumptions about the particle physics model. 
Instead, DAMA may simply test the caustic ring halo model by looking for the predicted 
big flow plateau in the recoil energy spectrum and the predicted annual modulation of that 
plateau. 
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APPENDIX A: COORDINATE SYSTEMS 



Here we enumerate the different coordinate systems used in the main body of the paper 
and state how they are related to one another. We use four coordinate systems: galactic, 
solar, terrestrial and "wake". All four coordinate systems are right-handed. 

The galactic coordinates are (xg,Vg, z g) with x G in the local direction away from the 
galactic center, y G in the direction of galactic rotation, and z G in the direction of the South 
Galactic pole. 

The solar coordinates are (X, Y, Z) with Z perpendicular to the plane of the Earth's 
orbit ("the ecliptic"), X in the direction of the Earth relative to the Sun at the time of the 
autumnal equinox (near Sept. 23), and Y = Z x X. The corresponding polar coordinates 
are called (6, $). 

The terrestrial coordinates are (a;®, y®, z®), where £® is the Earth's rotation axis, x® is 
in the same direction as X, and y® = £® x a;®. The corresponding polar coordinates (#®, 0®) 
are related to right ascension a and declination 8 by #® = | — 5 and 0® = a. 

The relation between terrestrial and solar coordinates follows from the fact that the 
ecliptic is inclined at 7 = 23.44° relative to the Earth's equator. We have 
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The relation between solar and galactic coordinates follows from the right ascension and 
declination of the Galactic Center (GC) and Galactic North Pole (GNP) [41]: a G c = 
17h 45.6m , 5 GC = -28° 56.2' , a GNP = 12h 51m , S G np = 27° 07.7'. This yields 
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To describe the effect of Sun's gravity on a particular dark matter flow, it is convenient 
to use "wake" coordinates. These are cylindrical coordinates (p,(p,z), with z axis is the 
direction of the flow velocity v far upstream of the Sun. Unlike ordinary cylindrical coordi- 
nates, however, we let p range from —00 to +00, and from to ir. This avoids the sudden 
change in value of the dark matter particles when they cross the z-axis. 
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APPENDIX B: THE CAUSTIC RING HALO MODEL 



The velocities and densities of the first 40 cold dark matter flows predicted by the caustic 
ring halo model at the Earth's location are listed in Table I. The velocities are with respect 
to the non-rotating rest frame of the Galaxy, and are given in the galactic coordinate system 
defined in Appendix A. Table I is an update of the table published in ref. [27]. The main 
difference is that the fifth flow is much more prominent in the new version. 

The caustic ring model explicitly realizes the fact that all cold dark matter particles 
lie on a three-dimensional sheet in six- dimensional phase-space. As a corollary there is a 
discrete set of flows at any location in a galactic halo. The number of flows is location 
dependent. There are caustic surfaces wherever the local number of flows changes by two. 
At the caustics the density of dark matter is very large. It was shown that two types of 
caustics, called inner and outer, are necessarily present in every galactic halo [22]. The outer 
caustics are simple fold (A 2 ) catastrophes located on topological spheres surrounding the 
galaxy. The inner caustics are elliptic umbilic (D_ 4 ) catastrophes located on closed tubes 
lying in the galactic plane. The inner caustics are usually referred to as "caustic rings of 
dark matter" . One caustic ring is associated with each flow in and out of the Galaxy. 

Observational evidence for caustic rings has been found in the distribution of rises in the 
Galactic rotation curve and in the presence of a triangular feature in the IRAS map of the 
Galactic plane [28] . The list of flows in Table I results from a fit to these observations. We 
do not justify Table I here, but merely use its flows as examples. 

The flows are labeled by an integer n — 1, 2, 3, .. and a binary symbol ±. In the caustic 
ring halo model, the first eight flows at the Earth's location (n = 1, 2, 3, 4) have large velocity 
components in the ±x G and +y G directions, whereas the remaining flows (n > 5) have large 
components in the ±z G and +y G directions. The two n = 5 flows have the largest local 
densities, by far. The reason is that the n = 5 caustic ring of dark matter, whose cross- 
section is allegedly revealed in the aforementioned IRAS map, is very close to us. One of the 
n = 5 flows has of order three times as much local density as all the other flows combined. 
So, it was named the "big flow". The density of the other n = 5 flow is much smaller than 
that of the big flow, but much larger than that of any of the other flows. The density 

of the big flow is large because of our proximity to a cusp in the nearby caustic ring. It 
is sensitive to our distance to that cusp. Because that distance is poorly known, there is 
a large uncertainty on the density of the big flow. Nonetheless it is clear that the big flow 
dominates over all other flows combined. 

For n < 4, d+ = d~ because the model is reflection symmetric. For n > 5, d+ and d~ 
are different in general. For n = 5 through 9, different estimates are given for d+ and d~, 
but for n > 10 the differences between d^ and d~ are not deemed sufficiently important to 
differentiate them. At present it is not possible to say, for each n > 5, whether d^ is the 
density of the flow of velocity v' 1 ^ or v" 1 ^. In particular this is so for n = 5. The velocities 
of the n = 5 flows are 

y 5 * = (470 y G ± 100 x G )km/s (Bl) 

in galactic coordinates. However it is not known whether v$ or v£ is the velocity of the big 
flow. 
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For n < 4, no entries are given for v^g because not enough information is presently 
available to estimate those velocity components. However the v^q velocity components are 
expected to be relatively small, < 30% of the corresponding velocity magnitudes 



n± 
G 



In the rest frame of the Sun, the flow velocities are 



where 



v Q = (-9x G + 232y G - 7z G )km/sec (B3) 

is the velocity of the Sun with respect to the Galaxy. In Eq. (B3) we assume that the Sun 
moves at 220 km/sec in the direction y G of galactic rotation plus 16.5 km/sec in the direction 
of galactic longitude / = 53° and latitude b = 25° [40]. The velocity components of the flows 
in solar coordinates, obtained from Table I by applying Eqs. (B2) and (A2), are listed in 
Table II. For the sake of definiteness we have set v^c = when computing the n = 1..4 
entries in Table II. 
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TABLES 



TABLE I. Velocity vectors v n± and densities d^ of the first 40 flows in the caustic ring halo 
model, in galactic coordinates. The flow of velocity vector v n± has density d^ or d%. 
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TABLE II. Densities and velocity vectors of the first 40 flows in the caustic ring halo model, 
in solar coordinates. 
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FIGURES 




FIG. 1. Schematic illustration of the effect of solar gravity on a CDM flow. Downstream of the 
Sun there are three flows, whereas upstream there is only one flow. The surface separating the two 
regions is the location of a caustic called the 'skirt'. It is a cone with opening angle equal to the 
maximum scattering angle of particles in the flow. 
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FIG. 2. The scattering angle S versus impact parameter b s , for the solar model of ref. [30] 
and incoming flow speed vq = 265 km/s. The dotted line is for a point mass Sun. 
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FIG. 3. The impact parameters b{ (i = 1,2,3) versus polar angle, denned relative to the 
incoming flow direction vq, at the Earth's distance from the Sun, for the solar model of ref. [30] 
and incoming flow speed vq = 265 km/s. The dotted line is for a point mass Sun. When \6\ exceeds 
the maximum scattering angle, there is only one flow and hence only one value (b\) of the impact 
parameter. When \6\ is less than the maximum scattering angle, there are three flows and three 
values (61,62 and 63) of the impact parameter. 
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FIG. 4. Flow speeds measured on Earth if the initial flow velocity is v 5 ~ , as a function of time 
of year, since Jan. 1. We neglect here the Earth's rotation and the eccentricity of its orbit. Flows 
2 and 3 exist on Earth from approximately March 30 until approximately Dec. 15. Their speed 
modulation is very small. Flow 1 exists throughout the year. Its speed modulation is of order 7%. 
The speed on Earth of flow 1 is largest near November 10 and lowest near May 10. 
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FIG. 5. Same as Fig. 4 but with initial flow velocity equal to v 5+ . Flows 2 and 3 exist on Earth 
from approximately June 13 until approximately March 7. The speed of flow 1 is largest near Jan. 
20 and lowest near July 20. 
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FIG. 6. Same as Figs. 4 and 5 but with initial flow velocity equal to v 10+ , whose direction is 
close to the ecliptic. The Earth passes near the flow's spike caustic on approximately Dec. 20. 
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FIG. 7. Flow densities di (i = 1, 2, 3) on Earth for flow 10+, as a function of time of year. The 
flow speeds for this case are shown in Fig. 6. The Earth passes near the flow's spike caustic on 
approximately Dec. 20. If = the densities of flows 1 and 2 diverge at that moment. In Figs. 
6 and 7, we used = 2.6 km/s, as an example. The Earth passes through the skirt caustic on 
approximately September 15 and March 25. However the associated divergences of the density of 
flows 2 and 3 are too feeble to show up on the graph. 
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FIG. 8. The function T(E r ), defined in Eq. (4.5), for the caustic ring model if the big flow has 
velocity , at those times of year when the speed of the big flow relative to Earth is largest(May 
10) and smallest (Nov. 10). For all n, the density d^ was taken to be that of the flow with velocity 
v"T. To establish the scale on the horizontal axis, we set the WIMP mass m x = 100 GeV and 
target nucleus mass tum = ^.w-proton with A = 73. For different values of m x and A, rescale the 
horizontal axis by a factor proportional to m^./m^ where m r is the reduced mass. 
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FIG. 9. Same as Fig. 8 except that now the big flow has velocity iJ 5 " 1 " and, for all n, is taken 
to be the density of the flow with velocity v^^. 
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FIG. 10. Event rate for the DAMA experiment in the 2-6 keV ee range, as a function of time of 
year, for m x = 50 GeV (top curve) and m x = 120 GeV (bottom curve), relative to the event rate 
on March 21. The velocity of the big flow was taken to be , as in Fig. 8. Spin-independent 
WIMP-nucleus couplings were assumed and the form factors were set equal to one. 
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